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Summary 

Although stringers are used primarily as stiffen- 
ers, they also can make damaged structures fail-safe 
or damage tolerant. Assessment of the damage toler- 
t ance of structures weakened by cracks is aided by 

knowledge of stress- intensity factors. In this pa- 
per, the stress-intensity factor is determined for a 
cracked orthotropic sheet adhesively bonded to an 
orthotropic stringer. Since the stringer is modeled as 
a semi-infinite sheet, the solution is most appropriate 
! for a crack tip located near a stringer edge. Both ad- 

herends are treated as homogeneous, orthotropic me- 
dia which are representative of many fiber-reinforced 
composite materials. It is assumed that the ad- 
herends are in a state of plane stress and the adhesive 
* is in pure shear. By use of Green’s functions and the 

complex variable theory of orthotropic elasticity de- 
veloped by Lekhnitskii, a set of integral equations is 
obtained. The integral equations are replaced by an 
equivalent set of algebraic equations, which is solved 
to obtain the shear-stress distribution in the adhesive 
: layer. From these adhesive stresses, equations for the 

stress-intensity factors at both crack tips are found. 

A parametric study is also conducted to deter- 
[ mine the sensitivity of the system to material prop- 

erties and specimen configuration. The parameter 
having the greatest influence on the stress-intensity 
i factors is the distance from the near crack tip to the 

edge of the stringer. Unless the crack tip is very close 
to or under the stringer, the stress-intensity factor is 
approximately that of an unstiffened sheet. How- 
ever, as the crack propagates beneath the stringer, 

J the stress-intensity factor decreases significantly. In- 

creasing the stringer stiffness or the adhesive stiff- 
| ness also results in a decrease in the stress-intensity 

factor. 

Introduction 

Because of their high strength and light weight, 
composite materials are finding increasing applica- 
tion in aerospace structures. In general, the struc- 
tural configurations used in composite airplanes have 
been very similar to the sheet-stringer construction 
currently used in metal airplanes. In metal airplanes, 
stringers are effective in making damaged structures 
fail-safe or damage tolerant. For this reason, the 
interaction of a cracked sheet and a stringer is an 
important problem which has been investigated by 
many authors. 

Romualdi, Frasier, and Irwin (ref. 1) considered 
two problems: a crack centrally located under a riv- 
eted stringer and a crack centrally located between 
l two riveted stringers. In their formulation of the 

problem, only the two rivets nearest the crack were 


considered effective. Sanders (ref. 2) considered the 
problem of a crack positioned symmetrically beneath 
a continuously attached stringer. He employed a 
shear-lag line stringer and assumed the sheet to be in- 
extensible in the direction parallel to the crack to ob- 
tain the stringer stress-concentration factor and the 
crack-tip stress-intensity factor. Greif and Sanders 
(ref. 3) removed Sanders’ previous inextensionality 
assumption and symmetry requirements. They con- 
cluded that the stringer-induced reduction in the 
stress-intensity factor was a localized effect and that 
the stringer stress-concentration factor was largely 
insensitive to sheet stiffness in the direction parallel 
to the crack. Bloom and Sanders (ref. 4) reported 
an analysis of a discretely attached stringer which 
may be broken. Poe (ref. 5) examined the problem 
of a cracked sheet stiffened by evenly spaced, riv- 
eted longitudinal stringers. The crack was assumed 
to originate at a rivet location and to grow beneath 
intact stringers. An extension of that work (ref. 6) 
considered the influence of broken stringers. 

Arin (ref. 7) examined the effect of the debond- 
ing of an infinitely long line stringer on the stress- 
intensity factor at the crack tip. He assumed that 
the stringer was adhesively bonded to an isotropic 
sheet along a line perpendicular to the crack. He 
concluded that the stringer exerts little influence on 
the stress-intensity factors unless it is quite close to 
the crack tip and the debond length is less than twice 
the crack length. 

Norris (ref. 8), using a complex- variable approach, 
represented the finite-width stringer with an array of 
line stiffeners. The stringer was divided into strips, 
and each strip was represented by a line stiffener 
attached to the sheet at discrete points, with no 
coupling between adjacent line elements. Because the 
analysis did not model debonding of the adhesive, 
it predicted unrealistically high loads in the line 
stiffener nearest the crack tip. 

Experimental work has shown that as the crack 
tip approaches the stringer, debonding can start and 
propagate through the adhesive. Although Arin 
(ref. 7) included adhesive debonding in his analysis, 
he modeled the stringer as a line, and thus was 
unable to account for partial debonding across the 
width of the stringer. Anderson, Hsu, and McGee 
(ref. 9) modeled partial debonding in their analysis 
of a cracked sheet reinforced by a bonded doubler. 
However, they were using a finite-element analysis, 
which would be very cumbersome and expensive to 
use when varying parameters such as crack length 
and debond area. 

A realistic analysis is needed that examines the 
effect of growing debond area on crack propagation 
in adhesively bonded structures. The purpose of this 



paper is to formulate such an analysis and to use it 
to study orthotropic sheets reinforced by orthotropic 
stringers. To meet this objective, the problem has 
been divided into two parts. In the first part, the 
problem is formulated with the assumption of a linear 
elastic adhesive with no debonding, and the stringer 
is modeled as a semi-infinite sheet bonded to an 
infinite sheet containing a crack. In the second part, 
the effects of adhesive nonlinearity and debonding are 
included. Although both parts have been solved, for 
brevity only the first part of the problem is presented 
here. 
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Pj,q 3 


QjkRjk 


Nomenclature 


a 

+a 
— a 

Aj , Bj , (sj , D j 


b 


c jk 

D 


half-length of crack, m 
right crack tip, m 
left crack tip, m 

constant coefficients in layer j 
stress functions, N/rir* 

distance from edge of stringer 
to center of crack, rn 

complex constants (y, k = 1,2) 

domain of integration 


S jk 

T 

1 3 


U 3 ’ V 3 


(b) (6) 

u\ \v\ 


u c , v c 


normal (mode I) component of 
stress-intensity factor, Pa^/m 

number of collocation points; 
number of cells in the domain 

complex constants for or- 
thotropic sheet (j = 1,2), 

1/Pa 

complex kernels used in 
displacement functions for 
layers 1 and 2, respectively 
(j,k= 1,2), 1/Pa 

complex kernels used in 
integral equations (j, k = 1,2), 
1 /(Pa-m) 

complex kernels used in stress- 
intensity factors (j = 1, 2, 3, 

4) 

x- and ^-displacements of layer 

3, m 

displacements in layer 1 due to 
crack-face pressure, m 

displacements in layer 1 due to 
concentrated forces, m 


E x> E y Young’s moduli in the x- and 

y-direetions, respectively, Pa 

Fp(w, z, a) complex- variable function 

used in displacement-potential 
functions 


(r) (c) 

, v) 


ir s , u ,s 


displacements in layer 1 due to 
distributed body forces, in 

displacements in layer 2 due to 
a pair of concentrated forces, 
m 


F h {w,a) 

complex-variable function used 
in calculation of left crack-tip 
stress-intensity factor 

F(?( w , z, a) 

complex-variable function used 
in displacement functions of 
layer 2 

Fn{z,w,a) 

complex-variable function used 
in stress potentials 

Fj(w , a) 

complex- variable function used 
in calculation of right crack-tip 
stress-intensity factor 

/oii /()2 

stress functions, m'^/N 

(• iy 

orthotropic shear modulus, Pa 


shear modulus of adhesive 
layer, Pa 

h j 

thickness of layer j\ m 

i 

= V-I 


u 2> v 2 

displacements in layer 2 due to 
distributed body forces, with 
respect to rotated coordinate 
system, m 

W J 

complex variable (j = 1, 2, 3, 
4), m 


body forces acting on layer j 
(. j = 1, 2), N/m' 4 

X' 2 ,Y< 

layer 2 body forces in rotated 
coordinate system, N/m 1 * 

x,y 

Cartesian coordinates, m 

2/0 

coordinates of load point, m 

x ,y 

rotated coordinate system, m 

Z 3 

complex variable (j = 1, 2), m 

zo 

location of load point, m 

Ax n , A y„ 

incremental distances 

(n = 1,2,..., N), m 
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L'xy : ^ yx 

roots of the characteristic 
equation (j — 1 , 2 ) 

Poisson’s ratios 

pressure applied to crack face, 


displacement potentials due 
to a pair of load points in a 
half-plane, m 



Superscripts: 


Pa 

(j)A k ) 

layer number (j, k — 1 , 2 ) 

T X) T y 

shear stresses in adhesive 
layer, Pa 

A bar over a 

symbol indicates a complex conju- 


gate; thus, if z = 

x -f iy , then z ~ x — iy. A prime 


stress potentials for or- 

on a symbol indicates reference to the rotated coor- 

thotropic materials, Pa 

dinate system, except in the stress and displacement 
potentials, for which a prime indicates differentia- 



displacement potentials due 
to a pair of load points in the 
cracked sheet, m 

tion. 



Formulation of the Problem 

Consider the sheet-stringer configuration shown in figure 1 . (Hereinafter, the semi-infinite sheet will be 
referred to as a stringer.) The stringer and the sheet are bonded together by an adhesive layer of constant 
thickness h%. Here the bond is assumed to be intact. The crack faces are subjected to a uniform pressure <jq 
and the stress state at infinity is zero, as shown in part A of figure 2 . The problem is solved with the following 
assumptions: 

1. The sheet (layer 1 ), the stringer (layer 2), and the adhesive (layer 3) are homogeneous and linearly 
elastic. 

2. The sheet and the stringer are dissimilar, orthotropic materials with principal directions of orthotropy 
being oriented parallel and perpendicular to the crack in layer 1 . 

3. The thicknesses of the sheet and the stringer are small compared with the in-plane dimensions, so that 
both layers are considered to be in a state of plane stress. 

4. The surface shear transmitted through the adhesive is assumed to act as a body force on the adherends. 

5. The thickness of the adhesive is small compared with the thicknesses of the sheet and the stringer; thus, 
the adhesive layer is treated as a shear spring rather than as an elastic continuum. 

Assumption 5 leads to the following continuity of displacement equations: 


^3 

ui -u 2 = —t x 


vi -v 2 = ~Q^ r y 


( 1 ) 


where u\,v\ and u 2 ,v 2 are the x- and y - components of the in-plane displacement vectors in layers 1 and 2 , r x 
and r y are the components of the adhesive shear stress, and hs and G 3 are the adhesive thickness and shear 
modulus, respectively. From assumption 3, the two sets of body forces (forces per unit volume) that act on 
layers 1 and 2 can be written as follows: 


*1 = - 


Tx 

h\ 





( 2 ) 


Figure 2 shows how the problem is broken into its component parts. The displacements in the sheet, shown 
in part B of figure 2, and in the stringer, shown in part C of figure 2, are determined individually as follows. 


Displacements in the Sheet (Layer 1) 

The forces acting on the sheet (layer 1 ) are shown in figure 3. The cracked sheet is subjected to the body 
forces X\ and Y\ distributed over the domain D and to the crack-face pressure oq. The problem shown in 
part A of figure 3 is equivalent to the superposition of the problems shown in parts B and C of figure 3. 
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c, 


If u^\v^ and denote the displacements in the sheet due to the loadings shown in parts B and 

respectively, of figure 3, the total displacements in the sheet can be written as 


u x 




Vl = v\ ' + V 




( 3 ) 


If we follow Lekhnitskii (ref. 10) and ignore an arbitrary rigid-body displacement vector, the displacements 
and due to the crack-face pressure shown in part B of figure 3, can be found. The complete 
expressions for and are given in appendix A. 

Determination of the displacements u\ ' and v\ ' due to the distributed body forces shown in part C of 
figure 3 starts with the known displacements due to a concentrated force located at the point zq = xq + iy q. 
Since the body forces Xi and Y\ are continuous functions of zq in the domain D, the displacements due to the 
concentrated forces can be used as Green’s functions to determine and Thus, these displacements 

can be written as follows: 


,(D 


(z,y) 


ff\Qu( 

D 


x, y, xo, yo) -Yj (xq, i/o ) + Qvi (x, y, xq, Vo) U (x 0 , yo )] dx 0 dyo 


\ c \x,y)= j J {Q‘2l{x,y,xo,yo)Xi (xq, yo) + Q 22 (x, y, xo, yo) Y\ (xq, 2 /d)] dx 0 dy 0 

D 


where Q\uQ\2iQ2i, and Q 22 are derived in appendix A. 


( 4 ) 


Displacements in the Stringer (Layer 2) 

The only forces acting on the semi-infinite stringer (layer 2) are the distributed body forces X^and Y 2 shown 
in part C of figure 2. In a procedure similar to that used in reference 11, the displacements due to a single 
pair of body forces X 2 and V 2 , located at point z q = Xq + iy f {) (see fig. 4), are used to form the displacements 
u? 2 and v 2 due to the distributed loading. Since the body forces X 2 and Y 2 are continuous functions of z^ in 
the domain, the displacements due to the concentrated force can be used as Green’s functions for the total 
displacements u 2 and v 2 . Thus, these displacements can be written as follows: 


u' 2 (jr\ y') = j J [i?i 1 (x' , y', Xq, y' 0 ) X 2 (xQ,y' a ) + R l2 (x', y, x' n , y' 0 ) Y 2 ( x' 0 , t/,',)] dx' 0 dy' 0 j 

1) 

v 2 {x' ■. y')= j J [^21 (*', y 1 , x'q, 2/0 ) x ' 2 (x' 0 , 2/0) + r 22 {x', y 1 , x(,, y' 0 ) Y 2 (x'q, 2/0)] dx'o dy[) 


( 5 ) 


The expressions for Bp, R\2,R2\i and R ? 2 are derived in appendix B. 

The displacements 112 and ^ 2 , required for the loading shown in part O of figure 2, are calculated by shifting 
and rotating the displacements vi 2 and v 2 according to the following relations: 


U2{x,y) = ~» 2 (x',y') v 2 (x,y) = u' 2 (x,y') 

where 


x' = y y' = b- x 

X' - Y Z Y’ = - X2 

Complete details 011 t he derivation of the displacements for the stringer are presented in appendix B. 


( 6 ) 


Integral Equations for r x and r y 

Substituting equations (3) and (6) into equations (1) yields the following system of integral equations for 
the unknown functions r x and r y \ 
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h 3 
Os' 


( x ) 2 /) + yy | ^y( x >y> x o,yo) + ^ ( x ',y', x o-yo) T x(x 0 ,y 0 ) 

D 

+ j^-{x,y,x 0 ,y 0 )- (x',y',x' 0 ,y' 0 ) r y (i 0 ,j/o)| dx o dy 0 = a 0 f 0l {x,y) 


h 3 , \,[[f \ / \ #12 / / / / / J 

— r y (x,y) + j J \ fov.zo.yo) - ( x ,y , x o.yo)' 


£> 


: ( x o,yo) 


+ ^ ( x ,y, x o,yo) + {z',y',x'o,y'o) ? y ( x o,yo)j dx 0 dy 0 = o- 0 /o2( x ,y) 


( 7 ) 


where 


/oi( x ,2/) = Real 




I 2 




- p<‘'4" (4 1 * - / 


- a* 


f 02 {x,y) = Real 


{ ^ 


(#) 


q\ x, t£ J I z\'> -\JzY r -a? 


,(1) 2 

'1 




x o = yo 

y'o = b - x 0 

x - y 

l i - =b- - x 


(4 1 ’ - i 


4 1)! - 


(8) 


The parenthetical numerical superscripts indicate layer number. 

As in reference 7, the kernels Qj and Rj^ (j,k = 1,2), which have logarithmic singularities, are known 
and are integrable in the domain D. Thus, equations (7) can be treated as Fredholm’s integral equations of the 
second kind and can be solved numerically by reducing them to a system of linear algebraic equations. If the 
integral signs are replaced by summations and collocation is used, the following system of linear equations in 
the unknowns r x (x, y) and r y {x,y) is obtained: 


N 


h 3 

G~ T x(xj,yj) + J2 [Sn ( Xj,yj,xo n ,yo n )T x (x 0n ,yo n j 

n= 1 

+ S 12 ( Xj,yj,x 0 n ,y 0 n )T y (zo„,yo„)] Ax n A y n = <r Q foi {■ Xj,Vj ) (j = l,N) 

hs N 

Q^ T y ( x jiVj) + Xj [ 5 2i {xj,yj,x 0n ,y 0n )T x ( x o„,yo n ) 

+ S 22 ( Xj,yj,x 0n ,y 0n ) T y (x 0n ,yo„)] Ax„A y„ = a 0 fo 2 ( xj,yj ) (j = 1,N) J 


( 9 ) 


where 


e Qn #22 

Sn = t; + 1^2 


S\2 = 


Qn _ #21 

h\ h-2 


(10a) 

(10b) 
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(10c) 


Q21 R \2 

S21 = TT " 

e Q22 Rn 

*22 = ~T t- - 7 — 

r»l /i2 


(lOd) 


and jV is the number of collocation points. The terms Q jk and R jk ( j,k = 1,2) are defined in appendixes A 
and B, respectively. 


The Stress-Intensity Factor 

In the present problem, because of symmetries in geometry and loading, the shear component of the stress- 
intensity factor is zero. The normal component is found by combining the effects of the crack-face pressure 
and the distributed body forces Xi(x,y) and Y\{x,y ) acting on layer 1. The stress-intensity factor may be 
expressed in terms of the unknowns t x and T y as follows: 


(«o) = 'Aix/a + J J [ T l ( T >y, x O,yo)T x (-To , </<)) +T 2 (x,y,x {) ,yo)T y (x 0 ,J/o)] dx 0 dy 0 

D 


(11) 


where 


T\ (x, y,x 0 ,yo) = cff [f, (44) + F, 

T 2 (x,y,x 0 ,y Q ) = [F/ (wj'4) + F/ («4 

s/w 1 - a 2 -(w- a 0 ) 

r I ( u\ a) — 

w — ao 


a) 

a) 


+c<; ) [f / ( 4 i) 

+4i } [*7 («4° 


«)+F/(4 1) 

o) + F t (u^ 


a )] 

“)] 


and ao is replaced by +a for the right, crack tip and by —a for the left crack tip. Complete details of the 
derivation of t he equation for st ress- intensity factors are given in appendix C. 

The solution of equations (9) gives the shear-stress distribution in the adhesive layer rj^x^y) and r y (x,y). 
Using these adhesive shear stresses in equations (11), we can determine the stress-intensity factors at either 
crack tip. The numerical solution to the problem is discussed in the next section. 


Numerical Solution 

A key item in the analysis is the method of in- 
tegration used for the system of integral equations 
represented by equations (7). In theory, these equa- 
tions require integration over the infinite domain. 
However, because of the complicated nature of these 
equations, a closed-form integration is probably im- 
possible. Consequently, the system of Fredholm’s in- 
tegral equations is solved with standard numerical 
techniques. This is done by dividing the domain into 
cells, assuming the unknown functions r x and r y to 
be uniform in each cell, and then using a numerical 
scheme to evaluate integrals. In this way the sys- 
tem of integral equations is reduced to the system of 
algebraic equations given by equations (9). Two im- 
portant aspects of the numerical integration are the 
domain of integration and the logarithmic singulari- 
ties in the kernels. 

In order to carry out the numerical analysis, the 
si/e of the infinite domain must be restricted. This 
restriction is possible because the disturbance caused 
by the presence of the crack decays with distance 


from the crack. A convergence study has been con- 
ducted to ensure that restricting the domain does not 
appreciably affect the accuracy of the results. Thus, 
the size of the domain D is determined iteratively, 
starting with a small, coarse-inesh domain and in- 
creasing the domain extent and mesh refinement un- 
til no significant changes occur in the quantities of 
interest. 

The second aspect of the integral equations re- 
quiring special attention concerns the kernels Sj^ 
which appear in equations (9). These kernels con- 
tain logarithmic singularities which must be evalu- 
ated separately. This is done by performing an exact 
integration of the singular portion of the kernels over 
the segment Ax n . In effect, this integration applies 
the load over a finite length rather than at a single 
point, thereby removing the singularity. 

Discussion of Results 

To gain insight into the behavior of the stringer 
panel, the analysis was performed for numerous 
combinations of material and geometric parameters. 
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Three different materials were assumed for the sheet. 
Their properties are given in table I. Material 4 , also 
shown in table I, was used to model the stringer. 
These materials span a wide range of properties and 
are representative of frequently used graphite/epoxy 
laminates (ref. 12). 

Figure 5 shows the normalized stress- intensity 
factor as a function of 6/a, the ratio of the stringer 
edge distance to the crack half-length. Here, the 
crack half-length a is held constant and 6 is increased. 
In this and all subsequent figures, the stress-intensity 
factors have been normalized with o^yfa , the stress- 
intensity factor for a cracked, infinite sheet with a 
uniform pressure ctq applied to the crack face. As 
shown in this figure, once the right crack tip is no 
longer beneath the stringer, the reduction in k\ is 
minimal. Also, the effect of the stringer on the 
left crack tip is much less than that on the right 
crack tip. This is to be expected, since the left 
crack tip is farther from the edge of the stringer. 
For large values of 6/a, both the right and the left 
crack-tip stress intensities approach 1.00 (the value 
when no stringer is present), again as expected. 
Changing the material properties of the sheet does 
not greatly affect the value of k \ , although the lower 
the longitudinal stiffness of the sheet material, the 
greater the reduction in the stress intensity due to 
the presence of the stringer. 

Figure 6 presents data similar to those shown in 
figure 5 , but in a slightly different manner. In this 
case, 6 is held constant while a is increased. The 
stress-intensity factor is not affected by the presence 
of the stringer until the crack tip is very close to 
the edge of the stringer, and no significant reduction 
in k\ occurs until the right crack tip is beneath the 
stringer itself. As in the previous figure, the stress- 
intensity factor at the left crack tip is not significantly 
affected by the stringer, and changing the material of 
the sheet does not have a large effect on k\. 

In figures 7 and 8, the effects of varying the stiff- 
nesses of the stringer and the adhesive are examined. 
In these figures, material 1 is used for the sheet and 
material 4 for the stringer; two ratios of 6/a, 1.2 and 
0 . 5 , are used. In these figures, only results for the 
right crack tip are shown, since any effects on the 
left crack-tip stress-intensity factor are minimal. 

(2) 

The effect of varying the stringer stiffness Ey }%2 
is illustrated in figure 7 . The normalized stress- 
intensity factor for the right crack tip is plotted as 
a function of the ratio of the stringer stiffness to 
the sheet stiffness. This figure shows that as the 
stringer stiffness approaches zero, the normalized k\ 
approaches 1.00, the solution for a cracked sheet with 


no stringer. As the stringer stiffness increases, the 
curves approach a lower bound on fcp This indicates 
that past a certain point, no further reductions in k\ 
can be obtained by increasing the stringer stiffness. 
The value of this lower bound on k\ is a function of 
the adhesive stiffness and 6/a. Changing the stringer 
stiffness has much less effect on the stress-intensity 
factor for the case of 6/a = 1.2 than for 6/a = 0 . 5 . 

The effect of varying the adhesive shear stiffness 
G^/hs is illustrated in figure 8. The normalized 
stress-intensity factor is plotted as a function of the 
ratio of the adhesive shear stiffness to the shear 
stiffness of the sheet. As the adhesive shear stiffness 
approaches zero, the value of k\ tends toward 1.00, as 
expected. As the adhesive shear stiffness increases, 
ki decreases. For 6/a = 1.2 the decrease is small, 
since the effect of the stringer is negligible. However, 
for 6/a — 0.5, the decrease in k\ is substantial. 
Although these points lie off the scale shown in 
figure 8, for very large values of G 3 (on the order 
of 10 10 MPa), the stress-intensity factor does not 
approach zero (as might be expected) but rather 
some nonzero minimum value. This lower bound is 
a function of the stringer stiffness and 6/a. Figure 8 

also shows that for values of G$h\/G^y /13 greater 
than about 4 . 0 , increasing the adhesive shear stiffness 
has little effect on the stress-intensity factor. 

Concluding Remarks 

This report has presented the formulation of an 
analysis for a cracked orthotropic sheet reinforced 
with an adhesively bonded orthotropic semi-infinite 
sheet. This configuration was assumed to represent a 
bonded stringer panel when the crack is located close 
to the edge of the stringer. The adhesive bond was 
also assumed to be intact. A parametric study was 
conducted to determine the sensitivity of the system 
to material properties and specimen configuration. 
The parameter having the greatest influence on the 
stress- intensity factor was the distance from the near 
crack tip to the edge of the stringer. Unless the crack 
was very close to or under the stringer, the reduction 
in k\ was minimal, and the farther the crack was 
under the stringer, the greater the reduction in k\. 
Increasing the stiffness of the stringer or of the adhe- 
sive also resulted in a reduction in the stress- intensity 
factor. 


NASA Langley Research Center 
Hampton, VA 23665 
April 11, 1985 
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Appendix A 

Derivation of Displacement Functions for Infinite Cracked Orthotropic Sheet 

This appendix presents the derivation of the displacement functions for the infinite cracked orthotropic 
sheet (layer 1). If we follow Lekhnitskii (ref. 10) and ignore an arbitrary rigid-body displacement vector, the 
displacements due to the crack-face pressure ctq (shown in part B of figure 3) can be written as follows: 


(Al) 


where 




Real 


I 




pi 1 ^^ (4 !) - \Al 1} -4*) - pi 1 V/ 1 (4° - 4 4 1)2 - « 2 ) 


u ( 1 1) -4 ,) 

»!*’ = [’! 1> ^ n (4” - W 7 *) - 4V," (4° - 44"^) 


4^ = 2 + 

„(*) _ _L_ 

- E {k) [Pj u *y ) 


a = crack half-length 

(3 = 1-2) 


(*) 1 
<?; = 


i 


(*) ..(*) 


(fc) (A:) 
VyxPj 


(j — 1,2; (/c) = Layer number) 
(j = 1, 2; (/c) = Layer number) 


(M (fc) (Jt) (A:) 

The terms /ij , ///, , , /i.^ are the roots of the characteristic equation 


,(k) 4 


fA k ) 
( 


K a 




,(Jt) ^ 


) '‘“ )i + ip = 0 ( ,m w*’) > °) 


where E x and E y are the Young’s moduli, G X y is the shear modulus, and v xy and u yx are the Poisson’s ratios 
of the orthotropic sheet. 

To determine the displacements for the distributed loading shown in part C of figure 3, start with the 
displacements for a concentrated force located at the point zq = xq -f- iyo. The displacements for a concentrated 
force can be written as 


u c {x,y) = 2 Real (4^) + P^ <t>2 (4^)] 

v c (x,y) = 2 Real 4/Vi (4^) + 4^2 (4^)] 


(A2) 


where 




Pi 

+ 


) 


JO 

P'2 


A\Fd (wW'zW^+BtFn (4 1) ,4' ) - a ) 
C X F D (4 1) ,4 ,) - a ) + D\Fd (^ 1) ,4 l) -«)] 

+ .4, log (zj L) - + Ci log (4» - W ,< 1) ) 


.(l) 


m - (1) Mi ' (1) [a,F /; (w\ l) ,4 l) ,a) + DyFo (4°, 4°, a) 


« ( 4 ”) 

+ C\ F i) +/?,F D (wV , '4' ) -«)] 

+ #1 log ( 4 ^ - w.44 + F>\ log ( 4 " - u>4^) 
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F£)(w, z, a) = log ^ \fvP‘ — a 2 \/z 2 — a 2 + u>z — a 2 ^ — log ^\/z 2 — a 2 + z^ 


(l) , (1) 

w\' =x o + Mi yo 

(l) (l) 

w 3 - x 0 - Ml yo 

Ai = cJJ ) ^ 1 + C' 1 ( J ) y 1 

Cl = cj^Xi - d$Yx 


(1) , (1) ' 
wy - x 0 + n j yo 

(i) (i) 

w\ =x 0 - n 2 yo 

Bi = c|J } Xi + C^Vi 

Di = cg^i - C^Yx J 


(A3) 


u o) r 
r U) _ L 

(l + ^ + vkp 

> 

°n “ . / 

27TZ ( /Li 
w 0’) 

?-&) (mP-m^)(mS j) -m? 

+ ^2^2^ + u *y 


°12 “ / 
27TZ ( /Li 

r (j)_ l 



2tti 1 fi 

Ji) 

Ai)- M2 


I 


). 


(A4) 


Details on the derivation of q ^zj 1 ^ and $2 ( 2 2 1 ' ) can be f° un d ' n reference 13. 

Since the body forces X\ and Y\ are continuous functions of zq in the domain D, the displacements u c and 
v c may be used as Green’s functions to determine and The displacements ^ and vj’ ^ due to the 

body forces X\ and Y\ distributed over the domain can be written as follows: 


u[ c \x,y)= J J u c (x, y, x 0 , yo) dx 0 dyo 
D 

= j J [Qn {x, y, I 0 . yo) (xo, yo) + Q12 (x, y, xo, yo) Yi (x 0 , yo)] dx 0 dy 0 

D 

v^(x,y) = J J v c (x,y, xo.yo) dx 0 dy 0 

D 

~ j j EQ21 (x, y, xo, yo) Ai (xq, yo) + Q22 (x, y, xq, yo) Yi (xo, yo)] dx 0 dyo 


(A5) 


where Q11, Q12, Q21, and Q22 are found by substituting equations (A2) into equations (A5) and then combining 
terms and simplifying. These expressions are shown below. 
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= 2 Real I 4^ [«!'>„ (»«", 4". a) + 4'X K". 4".«) 


Ml M*2 


- h', 4°. K'.-'J".-) 

7 „(i)„(i) ( r {\) , r (D\ r . X1 

- 2 \i n-n- ' "4/4"' -« 2 + 4")] 

M i M2 L V 


(Equation continued on next page.) 
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Appendix B 

Derivation of Displacement Functions for Semi-Infinite Orthotropic Sheet 

This appendix presents the derivation of the displacement functions for the semi-infinite orthotropic sheet 
(layer 2). The only forces acting on the semi-infinite sheet are the distributed body forces X 2 and I 2 , as shown 
in part C of figure 2. In a procedure similar to that used in reference 11, displacements due to a pair of body 
forces X 2 and Y 2 , located at point z f 0 = x f Q -f iy f 0 and oriented as shown in figure 4, can be written as follows: 


where 


where 


u s (x',y',x' 0 ,y' Q ) = 2 Real p ( , 2 V t 2) (4' 2) ) + pf (4 2) )] 
(x\y',x' 0 ,y' 0 ) = 2 Real (ef 2) ) + (^f)] J 


4 2) (4 2) ) = A ? 1o s (4 2) - u; i 2) ) + c 2 i°g (4 2) - <4 2) ) 

(2) {(M 2 >-4 2) )h^(4 2) -4 2) ) 


+ 


(2) 

M‘2 




+ 


c 2 i og (4 2 )-4 2) )] + (4 2) -4 2) ) 

B 2 log [z[ i] - «4 2) ) + log (4 2) - 


,( 2 ) 


)]} 


^ (4 2) ) = B 2 l„ E (.f>-»< 2 ))+D 2 log(J 


(2) - «.f ) 


( 2 ) 

F\ 


2 ) 1 (2) { ( 4 2) - 4 2) ) [m log ( 4 2) - 4 2) ) 

, ^2 


X (/i -2 log (4 J) - u)^) + D -2 log (4 2) - u/, 2 *)] } 


(Bl) 


(B2) 


zf = 4) + 4f )y 0 


O' = 1.2) 


and C 


>(2) ,.(2) c (2) 


1 1 


12 


21 


(2) 

„( 2 ) 


= 4) + M ( i 2) ?/0 


' 3 " = ~ x 0 + l l \ 2) 4) 

A 2 = c[?X! 2 + C^Y'l 

c 2 = -r|'f x' + 


( 2 ) / , ( 2 ) , 

w 2 = X 0 + M-i 2/0 

(2) _ _/ , ,.(2) 

’4 


w 


-i' 0 + [i 2 y[) 

B 2 = + 4? *2 

D2 = -4^ Y 2+4? y 2 


( 2 ) 


and are defined in equations (A4). Details on the derivation of <4 ^(4 all< ^ 


^ 2 ^) can be found in reference 13. 


As before, if the body forces X 2 and Kj are continuous functions of in the domain, the displacements u 
and n‘ s , given by equations (Bl), may be used as Green’s functions for the displacements u 2 and v 2 due to the 
distributed forces. Thus, the displacements can be found as follows: 
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u 2 { x '>y') = j J u * { x ',y'< x O,y'o) dx' 0 dy' 0 

D 

=//[*' l (x\ y', x'q, y' 0 ) X 2 (x' Q , y' 0 ) + R n (x', y', x' 0 , y' 0 ) Y 2 (x' Q , y' 0 )] dx' 0 dy' 0 

D 

v 2 (*'. y') = J J yS (*'> 2 /> * 0 - Vo) dx' 0 dy' 0 

D 

= J J [#21 ( x ', y\ X Q, y'o) X ’ 2 (i'o, 2/o ) + R22 { x , 2 /, *0, y'o) Y 2 K, 3 / 0 )] dz 0 d 2 /o 


where i?n, ^21> and #22 are found by substituting equations (Bl) and (B2) into equations (B3) and 

then combining terms. The complete expressions for #n, R\2, #2l ? and #22 are given below. 

Equations (B3) give the displacements of the semi-infinite sheet for the loading and the coordinate system 
shown in figure Bl(a). However, the displacements U 2 and V 2 required for the solution are defined with respect 
to the coordinate system and are due to the forces X 2 and Y 2 shown in figure Bl(b). Thus, from figure Bl, 


U 2 {x, y) = -v 2 (x, 2 /) V 2 {x, y) = u 2 (x', y') (B4) 

where 

x f = y y = b — x 

x f 2 = V 2 y 2 = -x 2 

The displacements in the semi- infinite sheet (layer 2) are now found by combining equations (Bl) and (B3) 
and substituting into equations (B4). Note that the material properties of layer 2 are defined with respect to 
the x f ,y f coordinate system (fig. 4). 

Shown below are the complete expressions for the functions #n, # 12 , R 2 i , and #22 > which appear in 
equations (B3). Further details on the formulation of these expressions can be found in reference 13. 



+ (4 2) - 4 2) ) <4? ['°g (4 2) - 4 2) ) - lo s (4 2) - 4 2) )] }) 
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#12 (x',y',x[ ) ,^) = 2Real(p ( 1 2) cg ^ 

+p [ ? c 22 [ ] °g (4 2) -«4 2) 


.(2) 


^cS 2 ) j^log (z {2) - W {2) ^j + log (. ^ 

"* '*>)+ log (*<*>-»?>)[ 


( 2 ) 

W 3 


■y 


{(<4 2, -4 2) 


l )ci 2 2 >[l.g( 2 < 2 >-«,' 2 ') + l„g( 2 ®-«>f) 


+ 


( 2 ) 

Pi 

4 2) - 4 2) " u ' 

(m 2 2> - M^) <4 2) lo S (4^ “ 4^) + *°g ( 2 [ 2) - ^4 2) )] } 

^ { („< 2 > - 4, 2 >) <?j 2 > [log (4 2 > - 4 2 >) + log (4 2) - 4 2) ) 

+ (4 2) - m ( i 2) ) ['°g (4 2) - <4 2) ) + lo g (4 2) - ^i 2) )] }) 


#21 ( x \ y\ x 'oi Vo) = 2 Real [log ^[ 2) - w| 2) ) - log ( 2 f 2) - w^)] 

+ 4 2)(7 2 2) [ lo S (4 2) - w 2 ] ) - lo S (4 2) - w 4 2) ) ] 


+ 


+ 


+ 


+ 


n 


( 2 ) 


{ ( A p) - 4 2 4 cj? [log (4 2 ) - *«) - log (4 2) - 4 2) y 


( 2 ) 

^2 


(2) 

Ml 


(4 2 > - 4 2 >) 4f [io g (4 2 > - »f) - io B (4 2 > - ) 


4 2) 


(2) (2) 
M l - M-2 


{ (M 2) - 4 21 ) A 1 ? b H 2) - *! 2) ) - ^ (4 s ’ - «?’)] 


(4 2 > - 4 2 >) 4? [i„g (4 2 > - »f) - i„ e (4 2) - 4 2 ')) }) 


ff'22 (- 


'>2/'mW)) = 2 #eal (r/[ 2 V;| 2) [log (z[ 2) - + log ^| 2) - w-’^)] 

/o\ /*>\ r / \ / b>>\i 

+ 4 


ai 2) 4 2) [log (4‘ 

J 2 ) 


M2) 

(2) \ I 

[ z 2 

w 4 ) 


+ (2)J (2) { (^ ( 1 2) - 4 2) ) C n [>°8 («1 2) - * 

/i 2 

+ (4 2) - 4 2) ) 4? [ lo s (4 2) - 4 l) ) + ! °g (4 2) - * 4 
"% U*?’ -<•?’) 4? M4 

.w-gf)]}) 


+ 


j(2) _ 

+ (4 2) - m! 2) ) 4? [iog (- 


Mi - M 2 

4 2) - «4 2) ) + io s (4 i 


’l^) + log (• 

■?’)]} 

•^ 2) - U’j 2) ) + log (4 


.(2) 

1 


,(2) 

'2 
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Appendix C 

Derivation of Stress-Intensity Factor Equations 

This appendix presents the derivation of the stress-intensity factor equations. In the present problem, 
because of symmetries in geometry and loading, the shear component of the stress-intensity factor is zero. 
The normal component is found by combining the effects of the crack-face pressure <jq and the distributed 
body forces X\(x , y) and Y\(x, y) acting on layer 1. For an orthotropic material, the normal component of the 
stress-intensity factor may be written as follows (ref. 14): 


= 2 V 2 — — (Tp - (1™ 0i (4^) + j J 0i ( 2 i^) dx odyo 

Z x ->a [ D 


(Cl) 


The function <j>\ , the stress potential due to the crack-face pressure ctq, is given by the following: 


(*!”) 


(1) 

a o 


do 


2 (mS 1) -4 1) ) v ^^i T ^ a2 


(C2) 


The function ^ the stress potential (derived in ref. 13) due to the body forces X\ and Yj, is given by 

the following: 


*< H") 1 2 V' F " 

Ml ~ M2 

+ B\Fn (« ( 1 ,) ,4 1) ,o)+C,F/ / ( 2 i 1) ^4 1) ,a) 

+ D\Fu 


+ C ' 


4 l > 


(C 3 ) 


where 


F H (z,w,a) = 


\f? — a 2 — \! u> 2 — a 2 — (z — w) 


(z - w) 

and A\, B\, C\, and Di are defined in equations (A3). Substituting equations (C2) and (C3) into (Cl), taking 
the limit, and combining terms yield the following stress-intensity factor: 

k\ = oo\fa - //*<’■ y, x 0 , yo) X\ (i 0 . yo) + T i (3;, y, *0, yo) U (^o. yo)] dx 0 dy» ( C4 ) 


where 


(1) N 
v) , a 

) +F] 

Key 

+ (?>!> 

(l) ' 

v\ , a 

)+F, 


j+41’ 

Fj(w , 

a) = 

\/ w 2 — a 2 — w; + a 

w — a 


Through the use of equations (2), equation (CM) may be written in terms of the unknowns r x and r y as follows: 
k\ =<T, )% /«+ J j [T 1 ! (*,y,x 0 , yo) r x ( 3 : 0 , yo) + T2 {x, y,XQ,yo)Ty (lo.yo)] dx n dyu (C 3 ) 
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Since the loading on layer 1 (as shown in fig. 3) is not symmetric about the ?/-axis, equation (C5) is valid only 
as z\ approaches a, that is, at the right crack tip. To find the stress- intensity factor at the left crack tip, the 
function <p ± must be found for the forces X\ and Y\ located in the left half-plane ( zq = — xq + x yo) • 

With the new form of 4 >\ same procedure as above is used to find /cj at the left crack tip. The 

function 4 >i does not change, since it does not depend on the location of the concentrated forces X\ and 

Y\. The stress- intensity factor at the left crack tip is thus found to be 


fc iL=-« = ^0^' 


y/a h 


i Jf * (x, y, I 0 , yo) T X (xq, 2/0 ) + ?4 ( x , y, xq, yo) r y (i 0 , 2 / 0 )] dx 0 dy 0 (C6) 

D 


where 


T3{x,y,xo,y Q ) = \^F f + F f +F f 

Ta (x,y, xo,yo) = cj^ [/> + Fp [V> + F f 


Fp(w,a) 


v w* — a 2 — (w + a) 
w + a 
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Figure 5 . Normalized stress- intensity factor as function of edge distance, a = 10 mrn; G 3 = 460 MPa; 
h 2 /h x = 0.5; h$/h\ - 0.005. 


No stringer 



Figure 6 . Normalized stress-intensity factor as function of crack length, b — 10 min; (1 2 — 460 MPa; 
h-i/hx = 0.5; h>Jh\ = 0.005. 
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